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Abstract 

The amplified spontaneous emission spectrum of the light field in the non- 
lasing supermode of two coupled semiconductor lasers is analyzed using lin- 
earized Langevin equations. It is shown that the interference between the 
laser mode and the fluctuating light field in the non-lasing mode causes spa- 
tial holeburning. This effect introduces a phase sensitive coupling between 
the laser field and the fluctuations of the non-lasing mode. For high laser 
fields, this coupling splits the spectrum of the non-lasing mode into a triplet 
consisting of two relaxation oscillation sidebands which are in phase with the 
laser light and a center line at the lasing frequency with a phase shift of ±7r/2 
relative to the laser light. As the laser intensity is increased close to threshold, 
the spectrum shows a continuous transition from the single amplified spon- 
taneous emission line at the frequency of the non-lasing mode to the triplet 
structure. An analytical expression for this transition is derived and typical 
features are discussed. 
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I. INTRODUCTION 



The carrier density dynamics of semiconductor lasers couple not only to the lasing modes, 
but also modify the amplified spontaneous emission into the non-lasing modes. It is there- 
fore possible to gain interesting insights into holeburning effects by studying the light field 
fluctuations in the non-lasing modes • The most simple system which can exhibit hole- 
burning effects is a laser array in which each of the individual lasers can be considered to 
operate in a single mode which is coupled to the local carrier density of that individual laser. 
A rate equation description for this type of device was developed by Winful and Wang |5||| 
and has been discussed in a more general context by the authors [Q. In the following, this 
model will be referred to as the split density model since the characteristic feature of the 
model is the separation of the carrier density into distinct density pools associated with the 
individual lasers in the array. Originally, the split density model was intended to describe 
instabilities in semiconductor laser arrays. Indeed, the model for a two laser array exhibits 
a region of instability with a limit cycle |J . However, the model can also describe the stable 
anti-phase locking between neighboring lasers in an array. This type of anti-phase locking 
is commonly observed in vertical cavity surface emitting laser arrays [|7|-[T0|]. In particular, 
the presence of carrier diffusion leading to a stochastical exchange of carriers between the 
lasers was found to enhance the stability of the anti-phase locking considerably [|]]. In the 
case of stable anti-phase locked operation, the special features of the model are still present 
in the dynamics of the quantum fluctuations causing amplified spontaneous emission in 
the non-lasing modes of the array. Specifically, the spatial holeburning represented by the 
carrier density difference between individual lasers is part of the fluctuation dynamics, mod- 
ifying the amplified spontaneous emission spectrum in the non-lasing modes according to 
this coupling between the light field and the carrier density dynamics |pT| . 

Experimentally, the symmetric non-lasing mode and the anti-symmetric lasing mode can 
be distinguished in the far field pattern. In particular, the spectrum of the non-lasing mode 
can be obtained in the center of the double-lobed far field where the anti-symmetric lasing 
mode vanishes due to destructive interference. At this very point, the non-lasing symmetric 
mode has an intensity maximum as a result of constructive interference. Assuming the 
lasing mode to be much stronger than the non-lasing mode, it should be possible to locate 
the center of the far field at the point where between the two maxima at the edges the total 
intensity is at its (local) minimum. If the far-field is observed off-center, a superposition 
of the laser spectrum and the spectrum of the non-lasing mode should be obtained. In 
this case it should be possible to distinguish the narrow laser line from the broad spectrum 
of the non-lasing mode. Spectral features of the fluctuations in the lasing mode such as 
relaxation oscillation sidebands should not be significant until the lasing mode intensity is 
much higher than the intensity of the non-lasing mode. Interference between the laser line 
and the fluctuations in the non-lasing mode are not to be expected since constructive and 
destructive interference are equally likely and consequently average out. 

In this paper we derive a general expression for the spectrum of the non-lasing mode in 
a two laser array. In section O, we introduce the split density model for two lasers operating 
in the antisymmetric supermode. The Langevin equation for the light field in the non- 
lasing mode is derived and the effects of spatial holeburning on the fluctuation dynamics 
are discussed. In section ITTTl the Langevin equation is solved and the analytic expression 
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of the spectrum is derived. General features of the spectrum are discussed. In section |TV], 
the spectra of four typical cases are presented, illustrating the different spatial holeburning 
features which may arise in different types of semiconductor array devices. 



II. SPLIT DENSITY MODEL OF A TWO LASER ARRAY 

A. Rate equations 

Generally a semiconductor laser array can exhibit a large range of dynamical effects 
involving far more than just two light field modes and two spatially separated carrier systems. 
However, the region of stable anti-phase locked operation relevant for at least some of the 
experimental realizations J7|-fn| can be described realistically by a split density model limited 
to only two light field modes and two carrier densities. The consistency of such a description 
with a more detailed model based on partial differential equations has been established in 



T2[ . We will use a version of the model which includes both the possibility of carrier diffusion 
and of a difference in the cavity loss rate between the symmetric and the antisymmetric 
supermodes. The rate equations read 

—Ek = ~Ni(l - ia)E x -(R + iu)E 1 - (~ - i~)£ 2 (la) 

j t N x = | - 7 Aq - F(N X - N 2 ) - 2wE* 1 E 1 N 1 (lb) 

jE 2 = ~iV 2 (l - %a)E 2 - (« + iu)E 2 - (~ - i^)E 1 (lc) 

j t N 2 = | - 7 JV 2 _ T(N 2 - Ni) - 2wE* 2 E 2 N 2 , (Id) 

where E x , E 2 are the field amplitudes of the individual lasers, normalized so that the inten- 
sities E^Ei, E£E 2 correspond to photon numbers in the cavity and N x , N 2 are the carrier 
densities above transparency in the active media of laser one and two, respectively, each nor- 
malized to represent the actual number of carriers in the laser. The total carrier injection 
rate \i is split equally between the two carrier pools. The linear gain coefficient w is then 
equal to the spontaneous emission rate into the laser mode. The rate of carrier losses due to 
electron-hole recombinations by spontaneous emission into modes not confined to the cavity 
is given by 7. The properties of the individual laser cavities are given by the frequency at 
transparency u and the cavity loss rate R. The carrier density dependent frequency change 
is represented by the linewidth enhancement factor a. 

The coupling between the two lasers is described in terms of the parameters T, s and 
Q. The carrier diffusion rate V represents diffusive carrier exchange between the two lasers 
which tends to equilibrate the carrier densities in the two regions of the gain medium. T may 
be derived from the ambipolar diffusion constant of the carriers, and the distance r 

between the two lasers, resulting in 

r = 4vr 2 ^2. (2) 

Note that the consideration of carrier diffusion effects represents an extension of the orig- 
inal model of Winful and Wang @||. The dissipative optical coupling term s represents 
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additional optical losses in the region between the two lasers, where the two field modes 
interfere. Constructive interference [E\ has the same phase as E 2 ) thus leads to increased 
losses, while destructive interference (phase difference of 7r and, therefore, opposite sign) 
leads to reduced losses. This coupling mechanism is itself sufficient to explain anti-phase 
locking 0,0- The particular physical explanation associated with this term in a real laser 
device depends on the type of waveguiding mechanisms used to separate the two lasers. In 
gain guided twin stripe arrays, the additional losses are due to internal absorption in the 
region of negative gain between the lasers. In vertical cavity surface emitting laser arrays 
fabricated by destroying the Bragg mirrors between the lasers [[7], this term represents the 
increased light emission from the cavity in that region. The coherent coupling of the two 
localized field modes by optical diffraction is represented by Q, the frequency difference 
between the symmetric and the anti-symmetric supermodes in the cavity at N\ = N 2 . 

Since this model represents a perfectly symmetric array of two lasers, it is useful to 
express the dynamics in terms of the symmetric and the anti-symmetric supermode, E + = 
l/V2(Ei + E 2 ) and £_ = l/V2(E l - E 2 ) respectively. Also, the carrier densities can be 
expressed in terms of the total carrier density N = Ni + N 2 and the carrier density difference 
A = Ni — N 2 . The equations then read 

j t E+ = ^N(l - ia)E + - (k+ + iuj + )E + + ^ A(l - ia)E_ (3a) 

jE_ = ^N(l - ia)E^ - («_ + iu_)E_ + ^A(l - ia)E + (3b) 

Ajv = fj, - 1 N - w(E* + E + + E*_E_)N - w(E* + E_ + E*_E + )A (3c) 

^-A = -(7 + 2r)A - w(E* + E + + E*_E_)A - w(E* + E_ + E*_E + )N, (3d) 

where k± = R ± s/2 and u± = u ± Q/2. In this formulation, the rate equations of the two 
laser array appear as rate equations of a single two mode laser with the additional variable 
A describing the spatial holeburning. 



B. Fluctuation dynamics 

If the effects of spatial holeburning are neglected by omitting all terms which depend 
on A, the stable solution is given by E + = 0, since k> + > Therefore, only the anti- 
symmetric supermode will show laser activity. In the presence of a non-vanishing spatial 
holeburning variable A, the stability analysis becomes more involved, as will be discussed 
below. However, there is always a stationary solution with 

A = 

N = 2k_/w 
E, = 



E_ 



1 /! 
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4k_ 


2w 



I exp[— i(uj_ + otK_)t], (4) 
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where Iq is the total laser intensity in terms of the photon number in the anti-symmetric 
supermode. The linearized dynamics of small fluctuations around this stationary solution 
may be separated into the relaxation oscillations of the laser mode and the fluctuations in the 
non-lasing mode and the holeburning parameter. This independence of the linear dynamics 
in E + and A from fluctuations in E_ and N is a consequence of the fact that E_ = A = is 
stationary for all values of E + and N. Only the relaxation and oscillation rates of the small 
fluctuations in E + and A are altered by E + and N. Therefore, the fluctuations in E + and 
iV appear only in the quadratic terms of the fluctuation dynamics, as products with A and 
E + , respectively. 

While the linearized fluctuations in the lasing mode describe the unaltered relaxation 
oscillations of a single mode laser, the fluctuations in the non-lasing mode are described by 

^-E+ = -(«+ - k_ + iu + + ian_)E + + -A(l - ia)sJT e -*(«--N™-)* (5 a ) 

^-A = -(7 + 2r + w/ )A - 2K^J7 (e- i(w - +aK - )t E* + + e +i(w - +a " e - )t £+). (5b) 

Note that the fluctuation dynamics can be written using 8E + = E + and 5 A = A because 
the averages of the non-lasing field mode E + and the holeburning parameterA are zero. 

Both the amplitude v% and the phase — (uo- + an )t of the average laser field appear 
in the fluctuation dynamics, representing the hole burning effects of the interference term 
between the weak non-lasing mode and the laser field. Thus the fluctuation dynamics are not 
independent of the phase of E + . The phase sensitivity can be represented in the fluctuation 
dynamics by expressing the complex field amplitude E + in terms of the component f\\ in 
phase with the laser field and the component f± which is n/2 out of phase with the laser 
field, i.e. 

E+ = (f\\ - if a) exp[-i(w_ + ««_)*]. (6) 

The fluctuation dynamics can now be rewritten as a linear matrix equation in a real three 
dimensional space. Adding the quantum noise terms Q± and Q\\, one obtains the linearized 
Langevin equation of the fluctuations in the non-lasing mode, 





wv^Fo/2 
-7 — 2T — wIq 





(Q\\ \ 


& 


Q± 




\ J 



(7) 



The quantum noise terms include quantum fluctuations of the light field vacuum entering 
the cavity and the dipole fluctuations of the gain medium. Here we have neglected the shot 
noise terms acting on the carrier density difference A. This assumption can be well justified 
in the case of R > 7, where within an order of magnitude of the threshold current the total 
number of carriers in the gain medium is much larger than the number of photons in the 
cavity . The relative shot noise in the carrier system will therefore be much lower than the 
photonic shot noise in the light field system. Note that the presence of relaxation oscillations 
is an indicator that this assumption is valid. The limit of its validity is reached when the 
stimulated emission overdamps the relaxation oscillations. 

The quantitative expression for the optical quantum noise terms Q\\ and Q± may be 
derived by comparing the corresponding intensity noise with the expected shot noise for 
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the quantized processes. The vacuum noise entering the cavity thus corresponds to the 
quantization of photons emitted from the cavity and the dipole noise in the gain medium 
corresponds to the quantization of photon emission from the gain medium into the cavity. 
Note that reabsorption of photons due to incomplete inversion may add extra noise. In the 
following, however, we assume perfect inversion in the gain medium by using the minimum 
noise terms required by quantum mechanics, 

(Q\\(t)Q\\(t + r)) = (Q±(t)Q±(t + r)) = k-8(t). (8) 

The fluctuation dynamics described by equation (|7|) depends on two relaxation times 
and two oscillation frequencies. In the case of weak laser light (Iq ~ 0), the light field 
dynamics are governed by the frequency Q representing the frequency difference between 
the supermodes and the damping coefficient s representing the loss rate difference between 
the supermodes. The interference effects between the laser field and the light field fluctu- 
ations in the non-lasing mode couple the holeburning parameter A to the field dynamics 
of the non-lasing mode. The total coupling strength is given by the relaxation oscillation 
frequency v /2/t_w/ - The damping of the holeburning parameter A is given by the sum of 
the spontaneous emission rate 7, the Diffusion rate 2T and the rate of stimulated emission 
wl . 

To illustrate the physical processes involved in the coupling between the spatial holeburn- 
ing parameter A and the field fluctuations in the non-lasing mode, it is useful to examine the 
real time dynamics of the fluctuations in h and f± for large intensities Iq. If v /2k_wJ 3> fl, 
relaxation oscillations dominate the dynamics of the fluctuations. However, since the hole- 
burning is only caused by the field fluctuation /y which are in phase with the laser field, any 
fluctuations in f± will not cause oscillations. Instead, f± relaxes exponentially with 

A « /|| « 

f± « f e- (s+an)t . (9) 

The effect of the frequency difference Q is suppressed by the fast relaxation oscillations and 
would not have any influence on the dynamics if it were not for the phenomena expressed by 
the linewidth enhancement factor a. The a factor describes waveguiding properties induced 
by spatial holeburning. If a carrier density difference A exists between the two cavities, a 
induces a change in the relative phase between the local fields E\ and E 2 . In terms of the 
supermodes, this means that the components of E_ and E + which are out of phase by 71/ 2 
are coupled. Therefore, A causes a change in f± via a by coherently converting part of the 
laser light amplitude into the non-lasing mode. In the presence of a strong laser field, even 
the negligibly small holeburning caused as Q rotates f± into h is sufficient to diminish f± 
by destructive interference with the light coherently transferred from the laser mode. 

In the relaxation oscillations caused by the spatial holeburning effect of interference 
between the laser light and in-phase fluctuations in the non-lasing modes, the a factor 
produces an out-of-phase component of a times the in-phase fluctuation, i.e. f± = af\\. 
Also, some of the field in f± is rotated into f\\ by the frequency difference Q to cause 
an effective undamping of the oscillations. Note that this is the exact counterpart of the 
additional damping in the exponential decay of f±. The real time evolution of the relaxation 
oscillations reads 
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/I, » / cos(^^)e^ +2 ™>+— ^/ 2 
/l « af cos(^2^wl t)e~^ +2r+wlo+s - an ^ 2 

A « -2^/^/o sin(^2Z^;t)e-^ +2r+w/o+s - aQ )*/ 2 . (10) 

If the dynamics are dominated by the relaxation oscillations, the exponentially damped 
dynamics of f± should produce a Lorentzian line at the laser frequency with a width of 
2s + 2aQ in the spectrum of the non-lasing mode, while the relaxation oscillation dynamics 
should form two sidebands shifted from the laser line by the relaxation oscillation frequency 
and having a width of 7 + 2T + w/ + s — aQ. If the a factor was zero, these sidebands 
would be fully symmetric. The main effect of a is to correlate the fluctuations in /y with 
the fluctuations in f± as expressed in equation (110). If the relaxation starts at f\\ = /o and 
f± — 0, the total relaxation dynamics are given by 

/||(t) = cos( v /2/ t _w/ t)e-^ +2r+w/0+s -^^/ 2 

f ± (t) = a (cos(V / 2^Iw7ot)e- (7+2r+w/o+s - an) ' /2 - e^ s+an ^ . (11) 

Figure p] shows this relaxation trajectory for a typical choice of laser parameters. The ques- 
tion of whether this trajectory contributes to a high frequency line or to a low frequency line 
depends on its curvature. The parts of the trajectory curved left rotate counterclockwise and 
contribute mainly to the low frequency line, while those parts curved right rotate clockwise 
and therefore contribute mainly to the high frequency line. Note that the sections of the tra- 
jectory contributing mainly to the low frequency sideband are significantly longer than the 
sections contributing mainly to the high frequency sideband. It is therefore to be expected 
that the high frequency sideband is lower than the low frequency sideband because of the 
a factor correlation between the relaxation oscillations along f± = ah and the exponential 
relaxation along /y =0. 

If afl is larger than 7 + 2T + s the relaxation oscillations may become undamped for 
intensities wIq < aVt — 7 — 2T — s. However, as the relaxation oscillations are slower at low 
intensities w/ , this effect does not necessarily occur at w/ = dVL — 7 — 2Y — s as predicted 
by equation (0). To find the correct stability boundaries, we will in the following examine 
the stability condition of the complete fluctuation dynamics given by S. 



C. Stability Analysis 



The stability boundaries associated with equation (p]) can be determined from the matrix 

-n 




-Wv^o/2 \ 
s — «wa/7q/2 
7 + 2T + w/ / 



(12) 



describing the linearized dynamics. At the stability boundaries, S must have at least one 
eigenvalue with a real part of zero. This condition is satisfied if there exists a frequency to 
with 
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det{S — iu} = 0. 



(13) 



The stability boundary for non-zero u corresponding to the undamping of relaxation oscil- 
lations by aVL can be found by separating the real and imaginary part of equation (|T3| ) and 
eliminating u. The condition for stability is then found to be 

an < s + 7 + 2T + w/ + — V f^ 2 + (s + 7 + 2r + w/ ) 2 ) . (14) 

ft_Wi v ' 

For the case of s = equation (|l4|) reduces to the expressions derived in previous studies 
@||. Since an instability is most likely to occur at low intensities Jo, however, a closer 
examination of the stability boundary is of interest in the present context. In particular, it 
is possible to determine the condition for no instability at any laser intensity, 



7 + 2T + s > an 1 1 + 2— - 2 J— (1 + — )(1 + — ) . (15) 
\ k_ \j k_ k_ cr J 

If the device parameters satisfy this condition, the anti-phase locked operation of the two 
laser array is stable over the whole operating range described by the two density model. 



III. SPECTRUM OF THE NON-LASING SUPERMODE 
A. General formula 

The Langevin equation (|7]) can be solved analytically in the frequency regime by de- 
termining the Greensfunction corresponding to the dynamical matrix S and applying it to 
the white noise input. The Greensfunction is obtained by inverting the matrix S + iu. In 
the two dimensional subspace describing only the optical field components f\\ and f±, the 
Greensfunction of the two density model is 



GM = 



1 1 s + iu n \ 

(s + iu) 2 + n 2 + M(an + s + iu) \~n~aM s + iu + M ) 

with M = (16) 

7 + 2r + w/ + iu v ' 

Note that M is a complex function of u and Iq which includes all parameters associated 
with the carrier dynamics. For M = 0, the Greensfunction is that of a cavity mode with 
a frequency of — n relative to the laser line and a damping rate of s. M increases with 
increasing laser intensity Iq, describing the transition from the cavity resonance to relaxation 
oscillations. 

Since the noise input given by Q\\ and Q± is white noise of intensity k_ in both compo- 
nents, the resulting fluctuations are 

l </||M7j_M) (I f±(u) \ 2 ) I 2n K ' K ' 



s 



The spectrum is given by the average intensities of the field as a function of the frequency 
uj. It can therefore be determined from the fluctuations in h and f± using 

I+(u) = {E* + {u)E + {u)) 

= (I /|| H | 2 ) + (| f±(fv) \ 2 )+i(f^(u)f ± (u)) - i(fl(cu)M<»))' (18) 

The spectrum of the non-lasing symmetric supermode I+(ou) for any given set of cavity 
parameters Q, s and any complex holeburning function M is then given by the general 
formula 

J+M = £ (19) 

/ 2(s 2 + (u - n) 2 ) + (1 - ia)(s - i{cu - Q))M + (1 + ia)(s + i{u - Q))M* + (1 + a 2 )M*M \ 
X { | (s + zw) 2 + Q 2 + M(s + iw + afi) | 2 J ' 



B. General properties of the spectrum 

To understand the physical processes associated with the spectrum described by the 
general formula given in equation ([OP, it is helpful to discuss the two limiting cases given 
by M = and M — > oo. For M = 0, equation ([T9|) describes the amplified spontaneous 
emission spectrum expected without spatial holeburning at a total carrier density pinned to 
N = 2k^/w, 

7 +< w >-° = T G- + ( .l + n)» ) • (20) 

Here, the total photon number in the non-lasing cavity mode is given by «_/s. In this case, 
no phase relation exists between the non-lasing mode and the lasing mode. In the opposite 
limit, i.e. M — > oo, only the component of amplified spontaneous emission out of phase with 
the laser mode remains undamped by the high stimulated emission rate and the spectrum 
shows a single line at the laser frequency uo = 0, 

T / x K- f 1 + a 2 \ .„ . 

Z7T \{s + ail) 2 + UJ Z J 

For M — > cxd the total photon number in the non-lasing supermode corresponding to this 
spectrum is + a 2 )/2(s + afi). For a = 0, this intensity is one half of the intensity at 
M — 0, reflecting the suppression of out-of-phase emissions. The linewidth enhancement 
factor a changes this simple situation however. It enhances the noise effect by a factor of 
1 + a 2 . This is the same mechanism which enhances the phase noise in the laser light and 
therefore increases the linewidth of the laser line. Moreover, a introduces an additional 
damping effect of afl, which makes this line wider than the one for M = and reduces the 
total photon number accordingly. 

The typical spectra pertaining to values of the complex rate M between these two ex- 
tremes are triplet spectra. The case described by equations fllOD features not only the 
centerline corresponding to M —>■ oo, but also a pair of relaxation oscillation sidebands 
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corresponding to the in-phase component of amplified spontaneous emission. In the limit 
of fast relaxation oscillations assumed in equations (10), the sidebands are centered around 



the relaxation oscillation frequency ujr = ±^2k_wI . Therefore, the sidebands can be de- 
rived from equation ([19]) by approximating the spectrum for large ±u;.r in the vicinity of 
uj = lor. The result of this approximation, including the lowest order asymmetry between 
the sidebands, reads 



k_ / 1 + a 2 



k_ II + a 2 + 2(ft + a(s + 7 + 2T + wI ))/y/2K-wl \ 
+ 2^ V (7 + 2T + w/ + s - aft) 2 + 4{u + y/2^wl ) 2 ) 
k_ ( 1 + a 2 - 2(ft + ajs + 7 + 2T + w/ ))/V2k_wV \ 
2tt \ (7 + 2T + wJ + s - aft) 2 + 4(cj - ^2k_w1 ) 2 J ' ^ ^ 

This triplet spectrum can be characterized by the total photon numbers in the three lines, 

(1 + a 2 W_ 



centerline 



2(s + aft) 



(l + a 2 )K_ / 2fl + 2a(s + 7 + 2r + w/ ) \ 

swfe " 4( 7 + 2r + wJ + s - aft) { (1 + a 2 ) v ^Iw7o J ' 1 J 

The sidebands are clearly damped by the stimulated emission into the laser line, wIq, as well 
as by spontaneous recombinations and carrier diffusion expressed by 7 + 2T. The increase 
in the stimulated emission rate is responsible for the suppression of the sidebands with 
increasing laser intensity. Note that the undamping of relaxation oscillations discussed in 
section || also appears in the total intensity of the sidebands, possibly raising the sideband 
intensities far above the intensity of the centerline at low laser intensities. 

The asymmetry of the sidebands can be traced back to two separate contributions, one 
proportional to the frequency difference between the modes ft and one proportional to 
a(s + 7 + 2T + w/ ). In both cases, the low frequency sideband intensity is greater than 
the high frequency sideband intensity. The first term can be understood as a remnant of 
the cavity resonance of the non-lasing supermode in the absence of spatial holeburning. 
The low frequency line continuously evolves from the line at uj = —ft for M = while 
the high frequency line only emerges as the laser intensity is increased. However, the low 
frequency line is also stronger for ft = 0. This is the effect explained in section ||.B which 
is caused by the a factor. Within the approximation assumed here, the sideband intensity 
ratio corresponding to the parameters given in Figure [l] is about two to one. 

While the general considerations given above already cover the qualitative effects which 
can be expected in the amplified spontaneous emission spectrum, the quantitative changes 
in the spectrum as spatial holeburning effects get stronger with increasing laser intensity are 
best illustrated and discussed on spectra pertaining to typical examples of device parameters. 
In the following, we will therefore discuss four cases, with either strong centerlines (ft < s) 
or weak centerline (ft > s) and either low carrier diffusion (7 + 2r «_) or high carrier 
diffusion (7 + 2T = «_/2). 
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IV. EXAMPLES OF SPECTRA FOR TYPICAL DEVICE PARAMETERS 

A. Diffusion rates and array size 

In the following we will discuss spectra assuming diffusion rates of T « 5 GHz to T rs 100 
GHz. Based on a typical ambipolar diffusion constant -D<-[iff in semiconductor lasers of about 
10 cm 2 s _1 the diffusion rate T is related to the separation r between the lasers according to 
equation (|2|) by 

/ r \ 2 _40GHz 
\l/xmy r 

Thus for GaAs based semiconductor lasers, the diffusion rates discussed correspond to arrays 
with separations between r = 0.6/im and r = 3/im. 

In this context the separations between the lasers of an array of about r = 3/xm represent 
diffusion effects which are much weaker than the spatial hole burning effects. This situation 
is to be expected for most of the semiconductor laser arrays realized today. Note that the 
spectra for higher diffusion rates and larger separations between the lasers are qualitatively 
similar to the examples given by the diffusion rates 7 + 2T = 10 GHz and 7 + 2T = 15 GHz. 
Consequently, these results also apply to array structures with much larger separations 
between the lasers. Below a separation between the lasers of r = 3fim the effects of diffusion 
drastically increase. In the following we investigate two representative cases of such strong 
diffusion effects with 7 + 2T = 50 GHz and 7 + 2T = 200 GHz. This corresponds to laser 
separations of about r = 1.25/im and r = 0.65/im, respectively. These values were chosen 
in view of a future realization of highly integrated and/or short wavelength semiconductor 
laser arrays. The spectra for 7 + 21" = 10 GHz and 7 + 21" = 15 GHz may thus be considered 
typical for the type of arrays already realized in experiments while the spectra for 7+21" = 50 
GHz and 7 + 2T = 200 GHz represent possibilities accessible only if further integration can 
be achieved. 

B. Triplet structure 

A typical triplet structure can be observed in devices where the frequency difference Q 
is smaller than the dissipative optical coupling s. The intensity of the centerline is then 
comparable to the intensity of the single line at M = 0. The carrier relaxation and diffusion 
rate 7 + 2T will suppress the sidebands, so a low value for this rate produces the strongest 
possible sidebands. The evolution of the spectrum as a function of laser intensity for the 
device parameters k_ = 400 GHz, s = 3 GHz, Q = 1 GHz, 7 + 2T = 10 GHz and a = 3 is 
shown in Figure ^| and || The laser intensity is given in terms of the stimulated emission rate 
w/o which can be converted to units of threshold current by dividing wIq by the spontaneous 
emission rate 7. Since 7 is usually close to one Gigahertz, realistic values for w/ range from 
about zero to ten Gigahertz. 

Figure |2| shows the evolution of the triplet structure with increasing intensity. Equation 
([22l) describes this evolution very well, except for the immediate vicinity of the threshold. 
Note that the asymmetry of the sidebands is extremely strong. At wio = 1 GHz, the choice 
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(24) 



of parameters corresponding to figure [I], the intensity ratio between the sidebands is about 
two to one, as predicted by the approximate equation (|22]) . 

Figure [3] displays an enlarged part of figure close to threshold. In this region, the 
triplet has not yet been formed. Between w/ = and w/ = 0.03 GHz the maximum of the 
single spectral line rapidly shifts to lower frequencies, increasing in intensity in the process. 
Between w/o = 0.03 GHz and w/o = 0.15 GHz, the line widens and finally forms a plateau 
on the low frequency side. At about w/ = 0.25 GHz, the plateau on the low frequency side 
splits from the center line to form the lower sideband. The high frequency sideband emerges 
in a similar fashion between w/ = 0.1 GHz and wIq = 0.5 GHz. 



C. Sideband suppression 

Since the center line is only damped by the sum of rates s + aQ while the sidebands are 
additionally damped by carrier diffusion, the sidebands will be suppressed in the presence 
of fast carrier diffusion. In this case, only a single line which gets frequency locked to the 
laser frequency as the laser intensity increases remains in the spectrum. This effect can 
indeed be seen in the spectra by choosing the same parameters as for the triplet structure 
(section [IV B ), k_ = 400 GHz, s = 3 GHz, Q = 1 GHz and a = 3, except for an increased 



carrier recombination and diffusion rate of 7 + 2T = 200 GHz. 

Figure |] clearly shows the evolution of a single line. However, the peak frequency is not 
just reduced to zero as the amplified spontaneous emission is phase locked to the laser light 
by the increasingly strong spatial holeburning effects. Instead, the peak even shifts to larger 
negative frequencies just above threshold. Between w/o = 0.5 GHz and wIq = 1.0 GHz the 
peak frequency is about 2 GHz, twice as large as the frequency at wIq = 0. This shift to 
higher frequency differences between the laser light and the non-lasing mode is similar to 
the one observed in the case of low carrier diffusion discussed above. However, the shift is 
not as rapid in the high diffusion case and instead extends to much higher laser intensities 
w/o- In the vicinity of w/ = 0, an analytical expression may be derived for the change in 
peak position by considering the coupling terms between the holeburning parameter A and 
the field components f\\ and f± in the dynamical matrix S as weak perturbations. The peak 
position uj p is given by the imaginary part of the eigenvalues of S. Near w/o = it shifts at 
a rate of 



d(w/ ) (7 + 2r) 

Note that this equation generally applies to all parameter sets. However, in the case of 
low carrier diffusion it is only valid extremely close to threshold. In particular, the shift of 
d(u p )/d(wIo) = 120 for the low diffusion case in figure |2| extends only to w/ = 0.03 GHz. In 
the low diffusion case displayed in figure |3] the shift of d(u> p ) / 'c?(w/ ) = 6 is the dominant 
feature of the spectrum up to intensities of w/o = 0.4 GHz, and even at w/o = 2 GHz the 
peak of the line is still 1 GHz below the laser frequency. These features in the evolution of 
the spectrum demonstrate how carrier diffusion suppresses the spatial holeburning effects. 
We note that a quantitative expression for the relative importance of spatial holeburning is 
given by the complex rate M introduced in equation (p!6|). It is the denominator of M which 
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includes the carrier diffusion rate 2T. Thus increased carrier diffusion reduces the absolute 
value of M, thereby suppressing spatial holeburning effects. 



D. Undamped sidebands 



If aVt is larger than 5 + 7 + 21", the undamping effect may cause an instability as discussed 
in section |IJ. If the instability is only narrowly avoided by choosing parameters close to the 
stability condition (|T5|), the sidebands should grow large and become narrow as the laser 
intensity passes the region of instability in the phase diagram. Figure [| shows the phase 
diagram for variable damping of the holeburning parameter, 7 + 2I\ and k_ = 100 GHz, 
s = 5 GHz, Q — 10 GHz, and a = 3. The horizontal line at 7 + 2T = 15 GHz shows the 
section of this phase diagram corresponding to the choice of parameters in figure |6|. The 
sidebands clearly have their maximum at a stimulated emission rate near w/o = 5 GHz, the 
point closest to the unstable region in the phase diagram. The linewidth of the sidebands 
in this region is close to 1 GHz, as compared to the linewidth given by 2s = 10 GHz. 
The maximal value of the total photon number in the non-lasing mode is roughly 240 at 
wIq = 5 GHz. For comparison, the total photon number at wio = is 20 and the centerline 
photon number for wio — > 00 is 14.3. The asymmetry of the lines is given by an intensity 
ratio of about three to one, e.g. at w/ = 5 GHz the photon numbers in the sidebands are 
approximately 180 in the low frequency sideband and 60 in the high frequency sideband. 



E. Sidebands in the presence of strong carrier diffusion effects 

If the carrier diffusion is increased, the undamping of the sidebands is suppressed. How- 
ever, the sidebands still remain the dominant feature of the spectrum if the line at wIq = 
is well separated from the laser line by Q > s. This effect can be illustrated using the same 
parameters as those in the case of undamped sidebands (section |1V 1J| ), k_ = 100 GHz, s = 5 
GHz, Q = 10 GHz and a = 3, except for an increased carrier damping and diffusion rate of 
7 + 2r = 50 GHz. As shown in figure [7|, the low frequency sideband continuously evolves 
from the amplified spontaneous emission line at threshold as intensity increases. Meanwhile, 
the weak phase locking effects cause a low intensity mirror image line to appear at the corre- 
sponding position on the high frequency side. Note that the frequency separation from the 
laser line increases from Q = 10 GHz at w/ = to almost 40 GHz at wio = 10 GHz. The 
relaxation oscillation frequency at wio = 10 GHz is about 45 GHz. This indicates that the 
line at wIq = 10 GHz is actually a relaxation oscillation line created by the carrier dynamics, 
while the suppression of the high frequency sideband is not only an effect of Q but also of 



the a factor as described by equation (|TTD. Indeed, the approximation given by equation 
(f22|) predicts an intensity ratio of twenty to one between the sidebands for the parameters 
used here. 



V. CONCLUSIONS 

We have demonstrated that the spatial holeburning caused by interference between the 
laser light in the antisymmetric supermode of a two laser array and the amplified sponta- 
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neous emission in the non-lasing symmetric supermode can give rise to relaxation oscillations 
which phase-lock the amplified spontaneous emission to the laser light. In the limit of high 
laser intensities, this effect splits the spectrum into a center line at the laser frequency 
and two sidebands. The center line represents the out-of-phase component of the ampli- 
fied spontaneous emission which does not interfere with the laser light. The two sidebands 
correspond to the relaxation oscillations between the in-phase component of amplified spon- 
taneous emission and the depth of spatial holeburning in the carrier distribution. 

The sidebands may become undamped by the linewidth enhancement factor a, which 
describes the conversion of laser light from the antisymmetric mode to the symmetric mode. 
This effect can be countered by carrier diffusion, which dampens the sidebands and may 
actually suppress them. Depending on the device parameters, it is therefore possible to 
alternatively find in the spectrum only the sidebands, only a single center line, or the full 
triplet. The continuous transition from a single line at the frequency of the symmetric 
supermode to any of the three possibilities can be described analytically using equation 

In all cases, the amplified spontaneous emission in the non-lasing supermode is strongly 
modified by the spatial holeburning effects. The spectra of the non-lasing modes above laser 
threshold are therefore quite different from the spectra at or below threshold described by 
linear optics. It should be possible to observe such spectra experimentally by measuring the 
spectrum of anti-phase locked laser arrays in the center of the farfield, where contributions 
from the anti-symmetric laser mode is close to zero. The type of spectrum observed will then 
allow a determination of the coupling strength and the type of coupling in the respective 
device. 
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FIGURES 



FIG. 1. Relaxation of f\\{t = 0) = 1, f±(t = 0) = for k_ = 400 GHz, s = 3 GHz, Q = 1 GHz, 
7 + 2r = 10 GHz, a = 3, and wIq = 1 GHz. The oscillations effectively correspond to a mainly 
counterclockwise rotation. 

FIG. 2. Amplified spontaneous emission spectra for k_ = 400 GHz, s = 3 GHz, 0=1 GHz, 
7 + 2r = 10 GHz and a = 3. (a) shows the contour plot of the spectrum as a function of laser 
intensity from wIq = to wIq = 2 GHz, while (b) shows the spectra at wIq = (no offset), 
w/ = 0.5 GHz (offset of 6/GHz), w/ = 1.0 GHz (offset of 12/GHz), and w/ = 1.5 GHz (offset 
of 18/GHz). 

FIG. 3. Amplified spontaneous emission spectra close to threshold for the same parameters as 
in figure 1. (a) shows the contour plot of the spectrum as a function of laser intensity from wIq = 
to w/o = 0.3 GHz, while (b) shows the spectra at w/o = (no offset), w/o = 0.1 GHz (offset of 
3/GHz), w/ = 0.2 GHz (offset of 6/GHz), and w/ = 0.3 GHz (offset of 9/GHz). 

FIG. 4. Amplified spontaneous emission spectrum for k_ = 400 GHz, s = 3 GHz, Q = 1 GHz, 
7 + 2r = 200 GHz and a = 3. (a) shows the contour plot of the spectrum as a function of laser 
intensity from wIq = to wIq = 2 GHz, while (b) shows the spectra at w/o = (no offset), 
w/ = 0.5 GHz (offset of 6/GHz), w/ = 1.0 GHz (offset of 12/GHz), and w/ = 1.5 GHz (offset 
of 18/GHz). 

FIG. 5. Stability boundary for k_ = 100 GHz, s = 5 GHz, Q = 10 GHz, a = 3 and variable 
carrier recombination and diffusion rates 7 + 2r. The diagonal line is the approximated boundary 
for sideband undamping given by 7 + 2r + w/o + s — a£l = 0. The horizontal line at 7 + 2r = 15 
GHz represents the choice of parameters in section 4 C and in figure 6. 

FIG. 6. Amplified spontaneous emission spectrum for k_ = 100 GHz, s = 5 GHz, 0, = 10 GHz, 
7 + 2r = 15 GHz and a = 3. (a) shows the contour plot of the spectrum as a function of laser 
intensity from w/o = top w/o = 10 GHz, while (b) shows the spectra at w/o = 3 GHz (no offset), 
w/ = 4 GHz (offset of 5/GHz), w/ = 5 GHz (offset of 10/GHz), w/ = 6 GHz (offset of 15/GHz), 
and w/ = 7 GHz (offset of 20/GHz). 

FIG. 7. Amplified spontaneous emission spectrum for k_ = 100 GHz, s = 5 GHz, = 10 
GHz, 7 + 2r = 50 GHz and a = 3. (a) shows the contour plot of the spectrum as a function 
of laser intensity from w/o = top w/o = 10 GHz, while (b) shows the spectra at w/o = (no 
offset), w/ = 2 GHz (offset of 0.2/GHz), w/ = 4 GHz (offset of 0.4/GHz), w/ = 6 GHz (offset 
of 0.6/GHz), w/ = 8 GHz (offset of 0.8/GHz), and w/ = 10 GHz (offset of 1.0/GHz). 
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